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Abstract. New developments in empirical analyses of the proton-proton differential cross section data at 
high energies are reported. Making use of an unconstrained model-independent parametrization for the 
scattering amplitude and two different fit procedures, all the experimental data in the center-of-mass energy 
interval 19.4 - 62.5 GeV are quite well described (optical point and data above the region of Coulomb- 
nuclear interference). The contributions from the real and imaginary parts of the amplitude beyond the 
forward direction are discussed and compared with the results from previous analyses and phenomenological 
models. Extracted overlap functions (impact parameter space) are outlined and a critical discussion on 
model-independent analyses and results are also presented. 

PACS. 13.85.Dz Elastic scattering - 13.85.-t Hadron-induced high-energy interactions 
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1 Introduction 

Elastic scattering is the simplest kinematic process in had- 
ronic collisions and the new data on proton-proton scat- 
tering to be obtained at the CERN-LHC by the TOTEIVQ 
Collaboration [1] has brought a renewed interest in the 
subject. In the theoretical (dynamical) context, high-energy 
elastic hadron scattering, as a soft process, constitutes a 
topical problem for QCD: perturbative techniques do not 
formally apply and nonperturbative approaches still de- 
pend on model assumptions, most of them unjustified [2] 
[3]. On the other hand, the experimental data available 
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can be well described by a wide variety of phenomenolog- 
ical models, but with different physical pictures [3| and 
therefore, an effective theoretical description of these pro- 
cesses, widely accepted, still remains an open problem. As 
recently stated, at this stage, "empirical parametrization, 
unbiased by any theoretical prejudice can be useful" [3] 
and that is the point we are interested in here. 

We have already treated model-independent analyses 
of elastic hadron scattering in a series of works jHISlinilZl 
I8ll9lll0). The goal is to extract from data empirical infor- 
mation on what is theoretically unknown, so as to justify 
some phenomenological inputs, getting insights for realis- 
tic model developments 6 and possible connections with 
QCD. 

The basic approach consists in the introduction of a 
model-independent parametrization for the scattering am- 
plitude, fits to the differential cross section data and em- 
pirical extraction of quantities of interest, as for exam- 
ple in an unitarized scheme, the overlap functions and 
the eikonal (impact parameter and momentum transfer 
spaces). However, that is not a trivial task for at least 
three main reasons: (1) the available data (differential 
cross sections) cover finite regions of the momentum trans- 
fer and the analytical integrations must be carried out 
through all the (physical) momentum transfer space; (2) 
experimental information on the phase of the amplitude is 
available only in the forward direction; (3) in a non-linear 
fit, the final result strongly depends on the choice of ini- 
tial values of the free parameters which are unknown. As a 
consequence, it is not possible to reach a unique solution, 
but only one or more possible solutions. 

The first problem can be treated by selecting only data 
with the largest values of momentum transfer available. 
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roughly above ~ 5 GeV^, and that limits the analysis to 
proton-proton (pp) elastic scattering in the energy region 
19.4 - 62.5 GeV, as we shall discuss. The other two prob- 
lems are crucial, demanding detailed analyses of different 
solutions, connected with different parametrization and fit 
procedures. 

We are presently investigating this subject and in this 
work we discuss new developments, with better statisti- 
cal results than those obtained in previous analyses 0[9]. 
However, we stress that we do not intend to give here the 
solution, but display and discuss what kind of solutions 
we can arrive with some specific methodology. As we shall 
discuss and explain along the text, only a global compar- 
ative investigation of different methods, parametrizations 
and solutions, extracting empirical information on what 
is common to all cases, can give strong support for theo- 
retical developments. We hope that this strategy and the 
results here presented can contribute with further devel- 
opments in the investigation of the so important inverse 
problems in high-energy elastic hadron scattering. 

The work is organized as follows. In Sect. 2 we re- 
call the main formulas connecting the physical quantities, 
scattering amplitude and the overlap functions (impact 
parameter space). In Sect. 3 we treat the improvements in 
previous analyses and in Sect. 4 we present the fit results. 
In Sect. 5 we discuss all the results, including comparisons 
with previous analyses and predictions from phenomeno- 
logical models, an outline of the emerging impact picture, 
our partial conclusions and some critical comments. The 
final conclusions and outlooks are the contents of Sect. 6. 



2 Physical quantities and overlap functions 

2.1 Physical quantities 

The physical quantities are expressed in terms of the com- 
plex elastic scattering amplitude A{s,q) — Re A{s,q) + 
i Im A(s, q), with — —t, where s and t are the Man- 
delstam variables. In this analysis we use as input the 
experimental data only on the differential cross section. 



the total cross section (Optical Theorem) 
<ytot{s) = 47r lmA{s,q = 0), 



(1) 



(2) 



and the p parameter 



ReA{s, q = 0) 



Imyl(s, q = 0)' 
From (1) to (3), the optical point is given by 



da 
d^ 



167r 



(4) 



2.2 Impact parameter space 

The impact parameter representation can be characterized 
by two basic functions, the Profile r{s,b) and Eikonal 
x(s,&), where b is the impact parameter. In the case of 
azimuthal symmetry we have [5] 



A{s,q) = i / bdbJo{qb)r{s,b), 
Jo 



(5) 



where Jq is the zero-order Bessel function and 

r(s,6) = (6) 

The Unitarity Principle, in the impact parameter space, 
can be expressed in terms of the total, elastic and inelastic 
overlap functions, 

Gtot{s,b)=Gel{s,b)+G,n{s,b), (7) 

which, in terms of the profile function reads [2] 



2Rer(s,6) = |r(s,6)|2 + G„(s,6), 



(8) 



where Gi„(s, b) can be interpreted as the probability of an 
inelastic event to take place at s and b. From Eqs. (6) and 
(8), 



G„(s,6) = l-e-2imxM)^ 



(9) 



where Im x(s, b) > and in the black disk limit G,;„ — ^ 1. 

As commented in our introduction, the model-independ- 
ent approach consists, essentially, in an empirical parame- 
trization for the amplitude A(s, q), fits to the da/dq^ data 
and the extraction of the overlap and eikonal functions [U 



3 Improvements in previous analyses 

In this Section we discuss the improvements introduced in 
the previous analyses by Avila-Menon [6 and Carvalho- 
Martini-Menon |9j, which concern three aspects: compi- 
lation and normalization of data at 19.4 GeV (Sect. 3.1), 
a novel unconstrained empirical parametrization for the 
scattering amplitude (Sect. 3.2) and new fit procedures 
and results (Sect. 3.3). 



3.1 Experimental data 

(3) 3.1.1 General aspects 



As already noted, the need of experimental information 
on the differential cross sections at intermediate and large 
values of the momentum transfer (up to and above ~ 5 
GeV^) limits the analysis to pp scattering at the highest 
energies with available data. Therefore, as in [S], we treat 
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Table 1. Compiled difTerential cross section data at 19.4 GeV, 
above the Coulomb-nuclear interference region: interval in mo- 
mentum transfer, number of points (N) and references. 



q-" interval (GeV^) 


N 


Reference 


0.0102 - 0.0315 


24 


Kuznetsov et al. [iSj 


0.021 - 0.66 


133 


Schiz et al. [H] 


0.075 - 3.25 


55 


Akerlof et al. [15] 


0.6125 - 3.90 


33 


Fidecaro et al. [16] 


0.95 - 8.15 


34 


Rubinstein et al. [17] 


5.0 - 11.9 


34 


Faissler et al. [18| 



here 6 sets of experimental data at the energies -y/s = 19.4 
(from the Fermilab and CERN-SPS), 23.5, 30.7, 44.7, 52.8, 
62.5 GeV (from the CERN-ISR). Each set includes the 
optical point, Eq. (4), and all the data above the region of 
Coulomb-nuclear interference, namely > 0.01 GeV^. 

Since the differential cross section data at large mo- 
mentum transfer, > 3.5 GeV^, do not depend on the 
energy at the ISR region, follovifing |6j[9] , we have included 
in each of the above five sets (ISR), the data at y/s — 27 A 
GeV (Fermilab), which cover the region 5.5 < < 14.2 
GeV^ (see |5] for a detailed discussion on this respect and 
a complete list of references) . 

The data from the ISR have been compiled, analyzed 
and normalized by Amaldi and Schubert and the en- 
semble of elastic pp data is the same used and quoted in 

except for the compilation and normalization of the 
data at 19.4 GeV, as discussed in what follows. 



3.1.2 Data at 19.4 GeV 

Previous analyses have shown that fit results at 19.4 GeV 
present some differences if compared with those obtained 
at the ISR region. These differences concern both the 
goodness of the fit and some extracted quantities. For ex- 
ample, evaluating the reduced chi-square, P^r degrees 
of freedom (DOF), and taking into account only the sta- 
tistical errors, typical results indicate x^/DOF ^ 2.76 at 
19.4 GeV and an average 1.48 ± 0.38 with the ISR ensem- 
ble [6|. Moreover, as shown in |6], the extracted eikonal in 
the momentum transfer space presents a zero (change of 
signal) with the ISR data, but not at 19.4 GeV. For those 
reasons, we have looked for all data available at 19.4 GeV 
and have investigated the necessity/possibility of a nor- 
malization. 

First, concerning the compilation, in addition to the 
data used in [5], and following [T^, we have included 
those obtained by Kuznetsov et al. [13] and Schiz et al. 
[14] ■ in the region of small momentum transfer. All the 
compiled data at this energy are listed in Table 1. The 
optical point, Eq. (4), is the same used in [B], namely 
77.66 ± 0.02 mbGeV-2, from atot = 38.98 ± 0.04 mb [H] 
and p = 0.019 ±0.016 

In a second step, without taking into account the sys- 
tematic errors, we have checked the consistence (normal- 
ization) of each data set with those at the nearby regions 
of momentum transfer (Table 1) and with the optical point 



(given above). In what follows we shall refer to each data 
set using the surname of the first author. Specifically, we 
have checked the consistence of the extrapolation of the 
data by Kuznetsov with the optical point, the data by 
Kuznetsov with those by Schiz and so on (Table 1), up 
to the data by Faissler. Figure 1 illustrates the sets at 
the diffraction peak (including the optical point) and in 
the largest region of momentum transfer. All sets are con- 
sistent from the optical point up to that by Rubinstein. 
However, as shown in the figure, the last set (Faissler) lies 
somewhat above the trend from the Rubinstein data. In 
this case we have taken into account the systematic error 
of 15 % in the Faissler data [18], scaling down the pub- 
lished results by a factor of 0.85. 



3.1.3 Data ensemble 

Summing up, in this work the ensemble of elastic pp data 
is the same used and quoted in [6] , except for the compila- 
tion and normalization of the data at 19.4 GeV. The new 
compilation at this energy includes the data by Kuznetsov 
and Schiz in the region of small momentum transfer (0.01- 
0.66 GeV^) and the normalization concerns a scale factor 
of 0.85 in the data by Faissler (largest region of momentum 
transfer). With this ensemble the first problem referred 
to in our introduction can be partially resolved, since the 
data cover the region up to ^ 12 GeV^ at 19.4 GeV and up 
to - 14 GeV^ at the ISR energy region (23.5-62.5 GeV). 



3.2 Unconstrained analytical parametrization 

In our previous work, we have used a constrained parametri- 
zation for the scattering amplitude, in the sense that the 
free parameters in the real part are also present in the 
imaginary part [B]: 



A{s,q) ^ 



with 



PCtot 



,(10) 



(11) 



E^tot 

where m < n and p, Utot are the experimental values at 
each energy. As shown in |B] , with this constraint the con- 
tribution from the real part of the amplitude presents a 
zero (change of signal) at small values of the momentum 
transfer, as predicted in a theorem by A. Martin for even- 
signature amplitudes |21| . 

In order to consider a more general and unconstrained 
form, we have introduced the following model-independent 
parametrization [7]: 



A {s,q) 



+ 1 ■ 



47r 



An 

n 

i=2 



E' 

4=2 



E«» 

i=2 



-diq^ 



' E c^^' 

J=2 



djq 



(12) 
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Fig. 1. Differential cross section data at 19.4 GeV at small 
(above) and large (below) values of the momentum transfer 
(see Table 1). 

where a^, bi (real part) and Cj, dj (imaginary part) are 
now independent free parameters. In the above formula, 
the parameters ai (real part) and ci (imaginary part) have 
been eliminated through Eqs. (2) and (3). As before, p and 
(Ttot are experimental inputs at each energy analyzed. 



3.3 Fit procedures 

Besides using the above unconstrained analytical function, 
our primary concern here is to develop fit procedures un- 
biased by any phenomenological (theoretical) assumption. 



In this section we discuss some critical aspects involved 
and after that, the methodology is developed. 

3.3.1 General aspects 

"Fitting nonlinear functions to data samples some- 
times seems to be more an art than a science." [22] 

Our purpose is to use the unconstrained parametriza- 
tion (12) to fit the differential cross section data, Eq. (1), 
at the six energies investigated. Since each data set com- 
prises ^ 150-300 points (covering 10 decades of experi- 
mental data) and the fits are non-linear in the exponen- 
tial parameters, we have no unique solution and we expect 
more than 10 free fit parameters at each energy. 

For comparison among different fit results we shall use 
the test for goodness of the fit. Since this test is based 
on the assumption of a Gaussian distribution of the data 
uncertainties |22i!22, we do not consider systematic er- 
rors of the data points, but only the statistical errors (ex- 
cept for the 0.85 normalization in the Faissler data at 19.4 
GeV). As a consequence, typical values of the per de- 
grees of freedom (DOF), for 100-200 DOF, are not so close 
to 1, as expected in terms of confidence intervals. However, 
even in this case, it represents a good test for compara- 
tive discussion on different results obtained with the same 
criterion. Anyway, if the systematic errors are included 
(for example in quadrature with the statistic errors), one 
obtains x^/DOF ~ 1, which seems to us unjustified on 
statistical grounds. 

Let us now discuss the other two main problems re- 
ferred to in our introduction, namely the non-linearity of 
the fit and the contributions from the real and imaginary 
parts of the amplitude, which are unknown beyond the 
forward direction. 

The non-linearity implies in a multimodal Chi-square 
function of the 2(m -I- n -1) free parameters, that is, a hy- 
persurface with multiple local minima and no direct access 
to a global minimum. As a consequence, any fit demands 
the choice of initial values of the parameters, which are 
also unknown and may lead to a local minimum without 
physical meaning and/or multiple solutions. 

Although parametrization (12) brings enclosed the full 
experimental information on the forward amplitude at 
each energy (namely Re A{s,0) and Im A{s^O) from atot 
and p data), the phase of the amplitude beyond this re- 
gion is experimentally unknown. Therefore, looking here 
for a completely unbiased fit, we shall not consider any 
constraint in initial values of the parameters associated 
with both the real and imaginary parts of the amplitude, 
nor put any limit in the number of parameters involved. 
Moreover, we are not interested to extract any energy de- 
pendence of the free parameter. The main point is to ob- 
tain a good fit on statistical grounds, at least better than 
those obtained in previous analyses 

The non-linear fits have been performed using the code 
CERN-Minuit [53], through successive runs of the MI- 
GRAD minimizer and with the confidence level for the 
uncertainties in the free parameter fixed at 70 % (vari- 
able up parameter depending on the number of degrees 
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of freedom). The error- matrix provides the variances and 
covariances associated with each free parameter. This in- 
formation is used in the evaluation of uncertainty regions 
in the differential cross section and all the extracted quan- 
tities, through standard error propagation procedures [22] . 

3.3.2 Methodology for initial values 

With all the above mentioned aspects in mind, we have 
considered two independent procedures for the choice of 
the initial values of the parameters, which we shall denote 
Methods 1 and 2. 

• Method 1 

In order to develop a first test on the response of our 
unconstrained parametrization to different initial values 
of the parameters, we based our choices on results already 
obtained in two other analysis, but excluding any con- 
straint or assumption involved, as explained in what fol- 
lows. 

The first one concerns the results by Avila and Menon 
[6], obtained with the constrained parametrization (11). 
For each data set we have used the values of the free pa- 
rameter (Table 5 in [B]) as initial values for the fits with 
our unconstrained parametrization (12), that is, the pa- 
rameters in the real and imaginary parts of the amplitude 
being now completely independent. 

The second one is related to a model-independent anal- 
ysis by Franga and Hama, in the same interval of energy 
[25] . The original parametrization is purely imaginary and 
the fits indicate a second dip in the region of large mo- 
mentum transfer. In our procedure, the values of the free 
parameters (Table 1 in [25]) have been used as initial val- 
ues in the imaginary part of our unconstrained amplitude. 
After a first run we have added one, two, three exponential 
terms in the real part until the fit reach the best x^/DOF 
(closest to 1), before starting to increase. 

We have also attempted analogous procedure with the 
parametrization by Amaldi and Schubert , but that 
form cannot describe the data at large values of the mo- 
mentum transfer (data at 27.4 GeV) and the quality of 
the fits was also very poor. 

Despite the differences in the structure of the two orig- 
inal parametrization and assumed conditions, the final 
results with the unconstrained parametrization and the 
above procedures were very similar (x^/DOF and the con- 
tributions from the real and imaginary parts of the am- 
plitudes). Therefore, we have selected for each energy the 
best statistical result (x^/DOF closest to 1) independently 
of the variant considered. The results are presented and 
commented in Sects. 4 and 5. 

• Method 2 

The second method is completely independent of the 
first one and has been developed in two steps. 

In the first step we have considered only the set at 52.8 
GeV, since it represents the best statistical information 
available (original data covering the region up to ~ 9.8 
GeV^ and with the addition of the data at 27.4 GeV, 
up to ^ 14 GeV^, totalling 245 points). The fit has been 
started with the data near the forward direction, <Q.l 



GeV^ and one exponential term in the imaginary part of 
the amplitude. In the logarithmic scale, the slope and the 
intercept of the straight line can be directly extracted and 
used as initial values for the fit. We have then extended the 
data up 0.5 GeV^, using the previous values as feed back 
and testing all possible number of exponential terms in 
the real and imaginary parts of the amplitude. Specifically, 
representing this number of terms as an ordered pair we 
have tested 

(real, imaginary) = (0, 1), (1, 1), (1, 2), (2, 1), (2, 2),..., 

until the fit result reach the smallest x^/DOF, before to 
increasing again. After that, with the same procedure, we 
have extended the data up to 1.5 GeV^ and then up to 
14.2 GeV^, using in each case the previous values as feed 
back and repeating the procedure. 

In a second step, once obtained the best fit at 52.8 
GeV, the values of the parameters have been used as initial 
values for the fits at the nearby energies, namely from 52.8 
to 62.5 GeV, from 52.8 to 44.7 GeV and then to 30.7, 23.5 
and 19.4 GeV. 

In all the above procedures, if the uncertainty in a 
given parameter was of the same order of, or greater than, 
its central value, this parameter was excluded and the fit 
performed again. We have also tested different confidence 
intervals for the uncertainties in the fit parameters, but 
the results were not as good as those with the fixed 70 % 
CL. Details on these results can be found in pB] . 



4 Fit Results 

The final fit results with Methods 1 and 2 are displayed in 
Tables 2 and 3, respectively. The visual description of the 
experimental data are quite similar in both cases and are 
illustrated in Figure 2, where it is also shown the uncer- 
tainty regions in the differential cross section, evaluated 
through standard error propagation from the fit parame- 
ters. 

The contributions to the differential cross sections from 
the real and imaginary parts of the amplitude, 

^ = 7r[RcAp, ^=^[ImAp, (13) 

are shown in Figures 3 (Method 1) and 4 (Method 2), 
where we display only the uncertainty regions and the 
experimental data. In the following we are going to discuss 
all these results. 



5 Discussion 

Table 4 summarizes the statistical results obtained here 
and those from previous analyses, for comparison. From 
that table, in terms of the x^/DOF (closest to 1) and 
taking into account all the six sets analyzed, we see that 
Method 1 presents the best statistical result. Except for 
the cases at 23.5 and 30.7 GeV the results with Method 2 
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Table 2. Fit results at each energy analyzed with Method 1: values of the free fit parameters and corresponding errors (aU in 
GeV"'^), number of points (A''), number of degrees of freedom (DOF) and reduced x^- 





-/s (GeV): 


19.4 


23.5 


30.7 


44.7 


52.8 


62.5 




(2l 


-0.03 


-0.05 


-0.022 


-0.035 


-0.043 


-0.022 






± 0.13 


d=0.40 


± 0.090 


± 0.094 


± 0.087 


± 0.097 




hi 


0.772 


0.93 


0.764 


0.858 


0.912 


0.774 


Re A 




± 0.034 


± 0.10 


± 0.046 


± 0.050 


± 0.034 


± 0.047 




0,2 


0.1825 


0.212 


0.3664 


0.5646 


0.7232 


0.8631 






± 0.0055 


± 0.023 


± 0.0047 


± 0.0063 


± 0.0053 


± 0.0046 




62 


1.760 


2.01 


2.644 


3.12 


3.249 


4.16 






± 0.061 


± 0.18 


± 0.092 


± 0.13 


± 0.099 


± 0.37 




Cl 


-34.08 


-32 


-157.30 


1244 


-128.8 


-140.1 






± 0.43 


± 16 


± 0.61 


± 22 


± 1.0 


± 2.7 




d^ 


2.5259 


2.63 


2 33972 


10.2433 


2 77889 


2.7678 






± 0.0024 


± 0.15 


± 00028 


± 0.0042 


± 00087 


± 0019 




C2 


392.00 


398 


3.802 


1460 


2.098 


1.87 






± 0.31 


± 11 


± 0.054 


± 10 


± 0.081 


± 0.19 




d2 


2.88166 


2.897 


8.51 


9.7419 


13.25 


14.4 






± 0.00033 


± 0.046 


± 0.18 


± 0.0034 


± 0.54 


± 1.3 




C3 


22.48 


9.4 


291.24 


-4.613 


-4.613 


226.3 






± 0.13 


± 3.1 


± 0.070 


± 0.070 


± 0.94 


± 2.6 




d-i 


4.023 


4.68 


2.43478 


2.283 


2.89792 


2.8746 


ImA 




± 0.013 


± 0.74 


± 0.00024 


± 0.012 


± 0.00051 


± 0.0012 




C4 


-373.86 


-369 


-219.53 


-2697 


-135.19 


-118.34 






± 0.27 


± 11 


± 0.19 


± 20 


± 0.28 


± 0.77 




di 


2.94863 


2.931 


2.68382 


9.9791 


3.2963 


3.2259 






± 0.00043 


± 0.045 


± 0.00035 


± 0.0020 


± 0.0015 


± 0.0043 




C5 


1.425 


1.48 


89.99 


1.57 X 10"^ 


56.30 


39.16 






± 0.075 


± 0.98 


± 0.13 


±0.84 X 10"^ 


± 0.18 


± 0.44 




ds 


11.35 


12.7 


2.9567 


0.399 


3.8306 


3.942 






± 0.56 


± 4.7 


± 0.0013 


± 0.052 


± 0.0052 


± 0.020 




C6 


2.98 X 10"^ 


1.7 X 10"^ 


1.15 X 10"^ 


5.702 


1.79 X 10"^ 


-1.22 X 10"^ 






± 0.59 X 10"^ 


±1.2 X 10"-^ 


±0.57 X 10"^ 


± 0.073 


±0.63 X 10"^ 


±0.72 X 10"^ 




de 


0.446 


0.405 


0.373 


2.435 


0.410 


0.377 






± 0.021 


± 0.073 


± 0.046 


± 0.013 


± 0.039 


± 0.056 




N 


314 


173 


212 


247 


245 


164 




DOF 


300 


159 


198 


233 


231 


150 




xVdof 


2.62 


1.09 


1.09 


1.99 


1.55 


1.18 



are also better than those obtained in the previous analy- 
ses. In this section we discuss the implication of the results 
with focus on the contributions to the differential cross 
sections from the real and imaginary parts of the ampli- 
tude (Sect. 5.1) and an outlining of the impact parameter 
picture associated with the overlap functions (Sect. 5.2). 
Based on this discussion and some other aspects we shall 
select Method 1 as our best result (Sect. 5.3) and then 
present some critical remarks (Sect. 5.4). 



5.1 Scattering amplitude 

The differential cross section data analyzed here are char- 
acterized by the well known diffractive pattern: the diffrac- 
tion peak, the dip region and the small monotonic decrease 
at large values of the momentum transfer. It is generally 
believed that this pattern is associated with a dominant 
imaginary amplitude, which is supposed to vanish only at 
the dip position, where the real part dominates filling up 



the dip. This standard picture seems to be inspired in the 
small value of the p parameter [q^ = 0) and in the fact 
that, at the energy region investigated, the p value goes 
through zero at the energies where the dip is sharpest 
P7] . This view seems also to be corroborated by Martin's 
theorem which states that the real part of the amplitude 
changes sign at small values of the momentum transfer 
pi] . suggesting an effective contribution only in the dip 
region. 

The above mentioned standard picture is in complete 
agreement with the results of our previous constrained 
analyses [Sl^ and also with the predictions from the great 
majority of the phcnomcnological models 3_ . However that 
is not exactly the case with the results presented here, as 
summarized in what follows. 

First we note that, from Tables 2 and 3, the numbers of 
exponential terms in Re A and Im A are rather different 
with Methods 1 and 2. In the former case Re A has 2 
terms and Im A has 6, and in the latter case, Re A has 
3 terms (except at 30.7 and 44.7 with 2) and Im A has 4 
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Table 3. Same legend as in Table 2 with Method 2. 





\/s (GeV): 


19.4 


23.5 


30.7 


44.7 


52.8 


62.5 




(Xl 


5168 


5489 


1730 


18.6 


51.7 


2780 






± 34 


±85 


± 210 


± 2.7 


± 6.2 


± 150 




bi 


7.43058 


7.4862 


7.1947 


8.09 


8.09 


8.0228 






± 0.00055 


± 0.0013 


± 0.0128 


± 0.24 


± 0.13 


± 0.0077 




(l2 


-1238.3 


-2182 


-1730 


-18.0 


-34.8 


-1226 


Re A 




± 93 


± 42 


± 210 


± 2.7 


± 3.4 


± 80 




&2 


7.1578 


7.3353 


7.2188 


6.25 


6.604 


7.787 






± 0.0019 


± 0.0027 


± 0.0031 


± 0.11 


± 0.072 


± 0.014 




13 


-3930 


-3307 






-16.2 


-1550 






± 33 


± 74 






± 5.2 


± 120 




b-i 


7.51291 


7.5767 






9.75 


8.190 






± 0.00068 


± 0.0020 






± 0.50 


± 0.012 




Cl 


-5200 


-570 


-1520 


-148 


-80 


-70 






± 120 


± 88 


± 170 


± 19 


± 10 


± 14 




rfl 


3.99471 


4.8389 


5.0777 


4.968 


4.909 


5.024 






±0.000072 


± 0.0037 


± 0.0012 


± 0.011 


± 0.019 


± 0.034 




C2 


5200 


579 


1530 


157 


90 


79 






±120 


± 88 


± 170 


± 19 


± 10 


± 14 




d2 


3.99892 


4.8645 


5.0891 


5.065 


5.077 


5.204 






± 0.00013 


± 0.0040 


± 0.0014 


± 0.012 


± 0.019 


± 0.031 




C3 


8.42 


0.0536 


0.0525 


0.0536 


0.0574 


0.0467 


ImA 




± 0.33 


± 0.0051 


± 0.0043 


± 0.0050 


± 0.0043 


± 0.0052 




d-i 


3.117 


1.026 


1.082 


1.062 


1.055 


1.008 






± 0.015 


± 0.056 


± 0.056 


± 0.056 


± 0.038 


± 0.057 




ca 


1.35 X 10"^ 


1.41 X 10"^ 


2.17 X 10"^ 


1.85 X 10"^ 


1.37 X 10"^ 


1.41 X 10-^ 






±0.075 X 10-2 


±0.66 X 10"^ 


±0.79 X 10"^ 


±0.77 X 10~^ 


±0.51 X 10"^ 


±0.64 X 10"^ 




di 


0.6027 


0.392 


0.431 


0.416 


0.388 


0.392 






± 0.0074 


± 0.046 


± 0.038 


± 0.042 


± 0.038 


± 0.044 




N 


314 


173 


212 


247 


245 


164 




DOF 


302 


161 


202 


237 


233 


152 




xVdof 


2.53 


1.77 


3.79 


1.87 


1.55 


1.18 



terms. From Figures 3 and 4, the uncertainty regions in 
the contributions from Re A and Im A are much larger 
with Method 1 than with Method 2. Taking into account 
these uncertainties the following features can be inferred 
in each case. 

- Method 1. 

Re A is maximum at =0, decreases monotonically 
and presents a zero (change of sign) in the region 1-2 
GeV^. 1mA has, in general, two zeros, one near the dip 
and the other around 4 GeV^, except at 62.5 GeV, where 
the second zero is not present. Im A dominates the peak 
and the region of large momentum transfer and Re A dom- 
inates the intermediate region, between the dip (~ 1.2 
GeV2) and _ 6 - 7 GeV^. 

- Method 2. 

Re A increases from = up to a maximum around 
0.3 GeV^ and then decreases to zero just above the dip 
position (without change of sign). Im A presents two ze- 
ros, the first around 0.5 GeV^ and the second at the dip 
position. Im A dominates the very forward peak and the 
region above the dip and Re A the peak. 

Indeed, our results with Methods 1 and 2 are rather dif- 
ferent from the above mentioned standard picture. Before 
discuss our conclusions (Sect. 5.3), let us outline some ex- 



Table 4. Values of the per degree of freedom (DOF) from 
previous analysis by Carvalho-Martini-Menon (CMM), Avila- 
Menon (AM), and those obtained in this work with Methods 
1 and 2. 



^ (GeV): 


19.4 


23.5 


30.7 


44.7 


52.8 


62.5 


CMM 19]: 


2.80 


1.20 


1.28 


2.13 


2.07 


1.51 


AM [B]: 


2.76 


1.20 


1.24 


2.05 


1.71 


1.22 


Method 1: 


2.62 


1.09 


1.09 


1.99 


1.55 


1.18 


Method 2: 


2.53 


1.77 


3.79 


1.87 


1.55 


1.18 



tracted results in the impact parameter space from Meth- 
ods 1 and 2. 



5.2 Impact picture 

By inverting Eq. (5), we extract the profile function and 
then the total, elastic and inelastic overlap functions, Eqs. 
(7) and (8), together with the corresponding uncertainty 
regions. Typical results obtained with Methods 1 and 2 
are shown in Figure 5, in the case of data at 52.8 GeV 
(largest set of experimental data). We see that despite all 
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1,5 



Fig. 2. Fit results and uncertainty regions from error propaga- 
tion: all region (above) and diffraction peak (below). Curves 
and data have been multiplied by factors of 10^^. 



the differences in the contributions from Re A and Im A 
the impact pictures are similar: the two results agree to 
within 20 %. 

Although not interested in the energy dependence of 
the free parameters, it will be useful to investigate the de- 
pendence of extracted quantities on the energy. However, 
it should be noted that the interval is relatively small, 
19.4-62.5 GeV and given the smooth variation of the phys- 
ical quantities in this interval (differential/total cross sec- 
tion and p) , the same smooth variation is expected on the 
extracted quantities. As illustration, we consider the ra- 



^ 1 1 1 1 1 1 


Method 1 












Black Disc Limit 
, — r— _ , 1 r : — 


in \ 


el \. 




b (fm) 

Fig. 5. Extracted overlap functions, Eqs. (7) and (8), at 52.8 
GeV, with the uncertainty regions from Method 1 and Method 
2. 



tio of the inelastic overlap function between two different 
energies. 



Rib) 



Gm(ysT, &)' 



(14) 



with = 62.5 GeV and y/si = 19.4 GeV (investigated 
interval). The results with Methods 1 and 2 are shown 
in Fig. 6, together with the uncertainty regions, by error 
propagation. Once more, despite the distinct contributions 
from Re A and Im A, we see that both methods predict 
a (standard) peripheral increase of Gin, although faster 
with Method 2 than with Method 1. 



5.3 Conclusion on the fit results 

As commented in the introduction, one of our primary 
concerns is to specify what kind of solution we can arrive 
with some specific methodology. We have considered two 
independent methods, based either in previous results and 
excluding constraints (Method 1) or by testing all kind of 
possible contributions in the real and imaginary parts of 
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• Method 1 
o Method 2 



" 



/ I + + + 

.fl' J. 





b(fni) 



Fig. 6. Ratio between the inelastic overlap functions, Eq. (14), 
at 62.5 and 19.4 GeV, with Methods 1 and 2. 



the amplitude and performing fits by regions of momen- 
tum transfer and energy (Method 2). From Table 4, except 
for the results at 23.5 and 30.7 GeV with Method 2, the 
reduced values are smaller and closer to 1 than those 
obtained in the previous analysis. Besides that, the con- 
tributions from Re A and Im A are distinct with Methods 
1 and 2 and both are rather different from the standard 
picture. Concerning this last point we make the following 
observations. 

- As explained in detail, we did not impose any theo- 
retical idea or constraint in the analytical parametrization 
or fit procedures either. 

- Up to our knowledge, the standard picture (dominant 
Im A with Re A filling up the dip) has never been proved 
or disproved. The arguments referred to in the beginning 
of Sect. 5.1 seem to us plausible, but not conclusive. 

- We have found that some extracted quantities of in- 
terest (exemplified here with Gin{s,b)) do not strongly 
differ with Method 1 or 2, and also from standard predic- 
tions of phenomenological models (we shall return to this 
important point in Sect. 6). 

Based on these observations, we understand that, al- 
though not standard, our empirical results suggest an- 
other picture as possible, in both statistical and physical 
grounds. However, we consider that Method 1 has led to 
the best results and that conclusion is based on the fol- 
lowing facts. 

(1) The reduced closest to 1 and smaller than that ob- 
tained in all our previous analyses and in all energies 
investigated (Table 4). To our knowledge that repre- 
sents the best global statistical result in the literature, 
if only statistical errors are considered. 

(2) Smallest uncertainty regions in the directly extracted 
quantities, which is connected with the correspond- 
ing variances, covariances and error propagation pro- 
cedures (Figure 5 is a typical example). 

(3) Re A presents a maximum ai — and its slope 
in this region is greater than that presented by Im A, 
which is consistent with a recent model-independent 



analysis by Kohara-Ferreira-Kodama on data at small 
momentum transfer j_l 2 1 . 

(4) The very poor statistical result obtained with Method 
2 at 23.5 and mainly 30.7 GeV, as compared with all 
the other results (Table 4). 

(5) The imaginary part of the profile function extracted 
with Method 2, at 30.7 GeV and at all other energies, 
has opposite sign: a reflection through the impact pa- 
rameter axis. 

Method 2, described in Sect. 3.3.2, seems to us a fully 
unbiased procedure and the above mentioned shortcom- 
ings in the final results were not expected. This intrinsic 
unbiased character of the method and our goal to dis- 
play what kind of solution can be obtained with a specific 
methodology were the reasons why we have included de- 
tails on the method and results in this work. Once selected 
Method 1 as our best result let us discuss some physical 
implications and critical aspects involved. 



5.4 Critical remarks 

What is novel in our selected result (Method 1) is the 
dominance of Re A at the intermediated values of the 
momentum transfer (2-7 GeV^) and also two changes 
of sign in Im A, except at 62.5 GeV (Figure 3). We note 
that the change of sign in Re A at small values of the 
momentum transfer is consistent with the prediction of 
the theorem by Martin [21], but that behavior was not 
imposed in the parametrization as done in fS]l2]. 

Despite the good fit quality with Method 1 (item (1) 
above) , some comments are in order regarding the numeri- 
cal values of the parameters at different energies, displayed 
in Table 2. First we note that the errors in the parameter 
fli are larger than the corresponding central values in all 
sets analyzed, which seems to violate the uncertainty crite- 
rion referred to in Sect. 3 (last paragraph). That, however, 
does not constitute an inconsistency since ai and also ci 
are not fit parameters: they have been eliminated through 
Eqs. (2) and (3) (note that the sum in parametrization 
(12) starts at i, j = 2) and have been included in Table 2 
only for completeness. This elimination also explains, in 
part, the apparently anomalous value of ci at 44.7 GeV, as 
compared with the values at all the other energies. How- 
ever, at 44.7 GeV, not only the value of ci but also of 
C2 and C5 seem anomalous and we understand this effect 
as associated with the lack of experimental data at the 
diffraction peak, in the region 0.3 < < 0.5 GeV^, which 
does not occur in the other sets (except at 52.8 GeV in 
the region 0.45-0.65 GeV^). That can also explain the high 
value obtained for the reduced at 44.7 GeV (compared 
with the other ISR sets), in all our analyses (present and 
previous), as can be seen in Table 4. 

The above mentioned and some others apparently anoma- 
lous values in different sets (for example C3 and C5) are 
also associated with our fit procedure and strategies: (1) 
method 1 was based in two different choices for the initial 
values of the free parameters (variants 1 and 2 referred 
to in Sect. 3.3.2) and the selection of the result was based 
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only on the reduced value, independently of the variant. 
That was possible because we did not identify any anoma- 
lous variation in the extracted quantities at the nearby 
energies; (2) with the Minuit code we did not impose any 
strategy in the search method, except for the starting val- 
ues and steps, and no limits have been imposed to any 
parameter, since we are looking for a completely unbiased 
procedure. 

On the other hand, it may be interesting to note that 
the above discussion applies only to the parameters asso- 
ciated with the imaginary part of the amplitude and not 
to those in the real part. In fact, from Table 2, all the 
values of the parameter bi are of the same order of magni- 
tude, lying approximately in the interval 0.76-0.93 GcV^^ 
and the parameters 02 and &2 present an increase with the 
energy in all the interval investigated. Despite these regu- 
larities, we stress that we do not intend here to investigate 
any energy dependence in the free parameters; our central 
point is to get the best result on statistical grounds. 

A distinct level of energy dependence concerns the 
extracted quantities and in that case, our results have, 
presently, limited applicability. The main reason is the 
relative small energy interval investigated (« 19-60 GeV), 
which also corresponds to the region just above the point 
where p changes sign and atot begins to increase, leading 
therefore, to several different possible uncontrolled extrap- 
olations. As a consequence of this reduced interval, the 
uncertainty regions, evaluated through error propagation 
from the fit parameters, generally overlap and do not allow 
detailed inferences on energy dependence at nearby sets. 
In this respect, the new data on the differential cross sec- 
tion from the LHC by the TOTEM collaboration, reaching 
the large momentum transfer region, will certainly shed 
light on this problem. 

On the other hand, even with these limited conditions, 
if we consider only the extreme energies of the interval in- 
vestigated, some trends can be inferred. As illustration we 
plot in Figure 7 the uncertainty regions associated with 
the real and imaginary parts of the profile function, at 
19.4 and 62.5 GeV, obtained with our selected Method 
1. Taking into account these uncertainties, we note a pe- 
ripheral increase in Re F and an increasing contribution 
in Im r through negative values, as the energy increases. 
In this interval both quantities present a change of cur- 
vature around 0.5 GeV, suggesting a change in the geo- 
metrical/optical structure, as predicted in the standard 
phenomenological context |:3:. 

We also note that if the new LHC data allow one to 
infer an energy dependence in our parameters, a predic- 
tive model-independent approach could be obtained, but 
that would demand the energy dependencies on Utot{s) 
and p{s), since they are input quantities in our analysis. 
However, that can be obtained, for example, from the de- 
tailed parametrization by the COMPETE Collaboration 

As commented in our introduction, some aspects of 
the inverse problem can only be treated through detailed 
analyses of different solutions, connected with different 
parametrization and fit procedures. However, some char- 




fa (fm) 



Fig. 7. Real and imaginary parts of the profile function at 19.4 
and 62.5 GeV, with Method 1. 



acteristics can already be extracted from our analyses, 
since they are all based on the inclusion of the data at 27.4 
GeV (large momentum transfer region) in the ISR sets, as 
explained in Sect. 3.1.1. As a consequence, all the empiri- 
cal results in these analyses indicate a smooth monotonic 
decrease of the differential cross section at high values of 
the momentum transfer, that is without any oscillatory 
or shoulder effect. Only to quote two typical examples, 
this result for pp scattering at 52.8 GeV favors the three 
pomeron model by Petrov and Prokudin ^29j , but not the 
hybrid model by Martynov and Nicolescu |30) . 

At last, some comments on the data sets presently 
available and those used here are in order. All our analyses 
have been based on the data sets published by experimen- 
tal collaborations; the complete list of references are given 
in Sect. 3.1.2 for the data at 19.4 GeV and in [B] for the ISR 
sets (Landolt-Bornstein Series) . Except for the normaliza- 
tion factor of the Faissler data at 19.4 GeV, we did not 
include the systematic errors in the analysis, only the sta- 
tistical ones, as has been justified in Sect. 3.3.1. Moreover, 
we did not make any kind of data selection, but included 
all the published results. For us, that represents a criterion 
that allows comparisons among different parametrization. 
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fit procedures and results, all based on the same data set. 
On the other hand, in the last years, some phenomenolo- 
gists have argued that the raw published data need selec- 
tions and also analyses on the systematic errors involved. 
Some compilations already exists, as the data set elabo- 
rated by Cudell, Lengyel and Martynov, discussed in [21] 
and available at "32^. We think it would be interesting to 
develop new data reductions with this data set and all our 
empirical parametrization. A comparative analysis on all 
the results may be useful (see comments on this respect 
in [30])- 

6 Final conclusions and outlooks 

The lack of a pure QCD description of the elastic scatter- 
ing data and the distinct physical pictures associated with 
the broad variety of phenomenological models, have been 
our primary motivation for exploring the inverse problem, 
as a source of empirical information, suitable for model 
developments and possible connections with QCD. In the 
context of an unitarized scheme (eikonal representation), 
results from previous analyses have already demonstrated 
the feasibility of this strategy, mainly related to empiri- 
cal information on the eikonal in the momentum transfer 
space [Blin]- Nevertheless, as discussed in detail, the main 
drawbacks concern the absence of a unique solution, which 
is a consequence of the non-linear fit and the fact that the 
contributions from the real and imaginary parts of the am- 
plitude, beyond the forward direction, do not constitute 
physical observables. 

Here our focus was in the development of a full un- 
constrained parametrization, fits and procedures, stressing 
the kind of solution that can be obtained with a specific 
methodology. In this context, despite the non standard 
picture related with the contributions from Re A and Im 
A (compared with phenomenological models and previous 
analyses), we have selected Method 1 as our best result. 
The novel feature is the dominance of the real part of 
the amplitude at intermediated values of the momentum 
transfer. 

At this point, the above mentioned drawbacks may 
suggest that we are faced with a messy situation, without 
conclusive directions/solutions and that the investigation 
of the inverse problem does not represent a way out. How- 
ever, the crucial point is to take these shortcomings as 
hints for developments and improvements: as commented 
in our introduction, we understand that an efficient strat- 
egy demands detailed investigation of different analyses 
and results. Since Re A and Im A above = do not 
constitute physical observables, the main point is to look 
for empirical information that do not depend on, or is 
not strongly connected with, these contributions. In that 
sense, the different results concerning our previous anal- 
yses (one zero in Re A and one zero in Im A) and those 
here selected (one zero in Re A and two in Im A), repre- 
sent only first steps in that direction. 

Summing up our strategies and outlooks are three fold: 

1. Global comparative analysis of model-independent pa- 
rametrization and fit procedures, with focus on good 



statistical results and distinguished by different con- 
tributions from Re A and Im A. 

2. Taking into account the uncertainty regions from er- 
ror propagation, selection and extraction of empirical 
properties that are common to all the analyses, namely 
that are not (or not so strongly) connected with spe- 
cific contributions from Re A and Im A. 

3. Use of these empirical information in the development 
of consistent phenomenological models, looking for pos- 
sible connections with QCD. 

At last we observe that, in addition to the previous 
analysis and that presented here, a third possibility 
has already been investigated through an almost model- 
independent representation for Martin's real part formula 
[33) (without the scaling property). This analysis has in- 
dicated two zeros in Re A and one in Im A [5, , what is dis- 
tinct from both our previous results and those presented 
here. Moreover, empirical information have been extracted 
which are consistent with the three analyses [33]. The in- 
vestigation along these lines are in progress and a global 
review on all our results are being prepared as a forthcom- 
ing work. 

We are thankful to J. Montanha for a critical reading of the 
manuscript. This research was supported by FAPESP (Con- 
tracts No. 09/50180-0, No. 07/01938-1, No. 07/05953-5) and 
CAPES, IFGW-UNICAMP. 
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Fig. 3. Method 1: uncertainty regions for the contributions to the differential cross sections from KeA(s, q) (grey) and lmA{s, q) 
(black) and the experimental data. 
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Fig. 4. Method 2: uncertainty regions for the contributions to the differential cross sections from KeA(s, q) (grey) and lmA{s, q) 
(black) and the experimental data. 



